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ABSTRACT 

The  use  of  optically  transparent  materials  to 
modulate  or  impede  laser  light  in  a  real  time  setting  is 
significant.  This  goal  stimulates  an  interest  in 
metalloporphyrins,  phthalocyanines,  naphthalocyanins 
and  their  derivatives.  The  most  desired  derivatives  would 
be  optically  transparent,  absorb  laser  light  and  also  be 
soluble  in  water.  Phthalocyanines  (PC)  with  their  stable, 
planar  square  structure  and  highly  delocalized  7i-electron 
system,  are  being  used  in  numerous  technological 
applications,  such  as  pigments,  in  chemical  sensors,  and 
more  recently  as  photosensitizers  for  photodynamic 
therapy.  The  nonlinear  optical  properties  (NLO)  of  these 
compounds  are  of  particular  importance.  Using  quantum 
chemistry  methods  we  study  both  visual  and  NLO 
properties  of  phthalocyanines  complexes  with  Si  as  a 
central  atom.  The  effect  of  hydrophilic  axial  substituents 
of  the  material  on  their  optical  properties  is  examined. 
We  are  mainly  interested  in  the  PEG  poly(ethylene 
oxide)  as  an  axial  substituent  because  it  is  used  in  our 
laboratory  to  synthesize  soluble  PC  complexes. 

1.  INTRODUCTION 

Phthalocyanines(PCs)  and  their  derivatives  have 
been  studied  extensively  since  the  beginning  of  the 
previous  century  being  important  dyes  and  pigments  due 
to  their  intense  blue  or  green  color.  In  recent  years,  PCs 
have  also  been  used  in  photocopiers  and  printers, 
photovoltaic  cells  gas  sensors,  nonlinear  optical 
modulating  devices,  photodynamic  therapy  agents  and 
many  other  applications  (The  Porphyrin  Handbook, 
2003).  The  PCs  display  interesting  properties  such  as 
high  thermal  and  chemical  stability,  efficient  light 
absorption  in  the  red  and  semiconducting  and 
photoconducting  abilities.  These  properties  are  result  of 
PCs  having  a  stable  macrocyclic,  conjugated  network  of 
7r-electrons  leading  to  high  electrical  polarizability  and 
rapid,  nonlinear  response  of  the  charge  density  to  the 
applied  intense  electromagnetic  fields.  The  large 
nonlinear  absorption  in  the  visible  spectrum  together 
with  the  ultrafast  response  time  and  easy  processability 
(Mckeown,  1998)  make  PCs  a  very  promising  modulator 
of  laser  beams.  Such  a  device  can  permit  the 


transmission  of  the  light  in  the  visible  region,  but  also 
strongly  reducing  high  intensity  laser  light.  As  such  the 
optical  modulation  properties  of  PCs  are  of  key 
importance  for  many  applications. 

Optical  modulation  with  PCs  was  first  reported  in 
1989  for  chloroaluminumphthalocyanine  (Coulter  et  al., 
1989)  and  since  then  many  other  PCs  with  these 
properties  have  been  identified  (Torre  et  al.,  2004) 
mainly  due  to  the  versatility  of  PCs  synthetic  organic 
chemistry.  The  new  PC  structures  can  be  made  by 
changing  the  chemistry  and  conformations  of  peripheral 
and  axial  substituents  as  well  as  by  varying  coordinating 
central  atoms  and  attaching  PC  rings  to  oligomeric  or 
polymeric  species  (Torre  et  al.,  2004;  Perry,  1997,  The 
Porphyrin  Handbook,  2003).  The  aromatic  7i-electron 
system  of  the  PC  can  be  expanded  resulting  in,  for 
example  napthalocyanies  (NCs)  that  also  exhibit  superior 
optical  modulation  properties  (Torre  et  al.,  2004). 

The  optical  modulation  properties  of  PCs  can  be 
explained  by  the  so-called  Reverse  Saturable  Absorption 
(RSA)  mechanism.  RSA  is  essentially  a  sequential  two- 
photon  absorption  leading  to  population  of  excited  triplet 
states  of  a  material  by  way  of  first  reaching  the  excited 
singlet  state.  The  optimal  RSA  material  would  have  an 
excited  triplet-triplet  cross  section  much  larger  than  the 
cross  section  of  the  ground  to  first  excited  singlet-singlet 
transitions.  It  is  also  necessary  that  the  intersystem 
crossing  from  the  lowest  singlet  excited  state  to  the 
lowest  triplet  excited  state  be  allowed  and  that  the 
excited  state  have  a  long  lifetime  (Perry,  1997).  A 
practical  optical  device  would  also  require  high 
concentration  of  soluble  RSA  material  in  the  optical 
beam.  This  in  turn  necessities  use  of  PCs  with  axial 
substituents  since  the  free  PCs  have  a  tendency  to  stack 
due  to  the  strong  intermolecular  interactions  between  the 
planar  PCs  rings.  Stacking  of  PCs  in  water  significantly 
affect  the  optical  spectra  and  degrades  their  optical 
modulating  capabilities.  Various  axial  substituents  have 
been  used  to  minimize  the  inter-molecular  stacking  of 
PCs  in  water  (Ngai  et  al.,  2001,  Shirk  et  al,  2000). 

In  this  paper  we  investigate  optical  absorption  as  a 
function  of  axial  substituents  in  PCs  with  Si  as  the 
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central  atom  (SiPC).  We  consider  hydrophilic 
substituents  such  as  poly(ethylene  oxide)  (PEG) 
oligomers  as  they  lead  to  a  material  soluble  in  water. 
Both  visible  absorption  (UV-vis)  and  non-linear  (RSA) 
are  studied  using  quantum  chemistry  methods. 

We  use  the  time  dependent  density  functional  theory 
(TDDFT)  (Gross  et  al.,  1996)  to  probe  the  excited  state 
properties.  This  method  has  been  previously  shown  to 
provide  reasonable  predictions  for  phthalocyanines  and 
metallo-phthalocyanines  structures  and  UV-vis 
absorption  (Day  et  al,  1988,  Ricciardi,  2001,  Liao  at  al., 
2004,  Wu  et  al.,  2004,  Liao  et  al.,  2005,  Zhang  et 
al.,2005,  Gunaratne  et  al.,2005).  Numerous  theoretical 
studies  have  also  been  performed  for  porphyrin 
complexes  (Baerends  et  al,  2002,  Rogers  et  al.,  2003). 

The  two  main  excitations  peaks  of  porphyrins,  the  Q 
band  in  visible  range  and  the  B  (or  Soret)  band  in  the  UV 
region  has  been  explained  by  the  “four-orbital  model” 
introduced  in  1963  by  Gouterman  (Goutermann  et  al., 
1963).  According  to  this  model,  the  Q  band  originates 
from  electron  excitations  out  of  the  HOMO  and  into  the 
LUMO  and  LUMO+1,  while  the  B  band  is  due  to 
transitions  from  HOMO-1  to  LUMO  and  LUMO+1. 
Later  studies  have  proposed  a  more  detailed  explanation 
involving  more  than  “four  orbitals”  because  the  quasi¬ 
degeneracy  of  orbitals  close  to  the  HOMO  is  removed. 
(Toyota  et  al.,1997).  In  the  case  of  metallo- 
pthalocyanines,  transition  metal  orbitals  are  also 
involved  in  the  excited  states  (Edwards  and  Gouterman, 
1970). 

We  are  not  aware  of  any  published  theoretical 
investigation  for  non-linear  RSA  properties  of 
phthalocyanines.  Unlike  study  of  UV-vis  spectra  which 
requires  calculation  of  singlet-singlet  excitations  (S-S), 
prediction  of  RSA  requires  calculation  of  triplet-triplet 
(T-T)  absorption  and  ultimately  excited  singlet-triplet 
cross-section.  A  few  investigations  have  been  done  on 
excited  triplet  states  for  porphyrin  derivatives  and  zinc- 
porphyrin  (Nguyen  et  al.,  2000)  as  well  as  meso- 
diaryloctaalkyl  porphyrins  (Kyrychenko,  2002). 

2.  COMPUTATIONAL  METHODS 

The  conformational  analyses  of  the  molecular 
systems  described  in  this  study,  including  structural  and 
orbital  arrangements  as  well  as  property  calculations, 
were  carried  out  using  a  variety  of  computational 
techniques  for  comparative  purposes,  using  GAMES S 
(Schmidt  et  al.,  1993),  GAUSSIAN03  (Frish  et  al.,  2004) 
and  DMol  (Delley,  2000).  Choice  of  optimal 
methodology  was  made  by  performing  computations  on 
the  lead  target  compound,  using  B3LYP,  PBE,  LDA  and 
PW91  functionals  in  combination  with  a  variety  of  basis 
sets  including  6-31G(nd,mp)  and  6-31  lG(nd,mp)  [n  =  1, 


2;  m  =  0,  1].  Molecular  structures  were  optimized  using 
the  DMol  program,  with  the  double  numerical  basis  set 
at  the  PBE  Hamiltonian  level.  We  used  the  delocalized 
internal  algorithm  for  geometry  optimization.  Although 
this  is  a  highly  efficient  algorithm  it  can  only  find  a  local 
minimum  close  to  the  initial  structure. 

Therefore  in  order  to  explore  the  larger 
configurational  space  of  the  rather  flexible  PEG 
oligomers  attached  to  PC,  we  have  used  molecular 
dynamics  (MD)  technique  to  generate  realistic,  initial 
structures  for  DFT  geometry  optimization.  The  initial 
structures  of  compounds  with  PEG  in  water  were 
obtained  from  MD  simulations  of  the  molecules 
immersed  in  water,  represented  by  a  unit  cell  of  500 
water  molecules.  The  structures  of  SiPC  were  kept 
frozen  at  their  DFT-optimized  geometries  while  PEG 
conformation  and  water  were  subjected  to  heating  the 
system  to  900K  for  50ps  followed  by  annealing  in  steps 
of  5ps/50K  and  finally  equilibration  at  300  K.  Discover 
program  with  the  COMPASS  forcefield  and  the  NVT 
Berendsen  thermostat  was  used.  The  annealing  procedure 
was  found  necessary  to  allow  for  proper  distribution  of 
water  molecules  in  the  unit  cell.  After  100  ps  of 
equilibration  the  snapshots  from  such  simulations  were 
taken  and,  after  further  DFT  geometry  optimization,  used 
in  calculations  of  excited  states.  We  have  identified 
several  conformers  within  the  energy  range  of  4  kcal/mol. 
The  lowest  and  the  highest  energy  conformers  were  used 
in  subsequent  calculations.  We  have  also  retained  some 
water  molecules  within  the  van  der  Waals  radius  of  the 
PC  structure.  This  approach  allows  us  to  generate  various 
possible  structures  of  the  PEG  polymer  attached  to  the  Si 
atom  of  PC.  The  optimized  structures  are  available  as 
Supplementary  material. 

The  UV-Vis  spectrum  was  identified  from  singlet 
excitations,  while  the  RSA  spectrum  was  found  from  the 
triplet-triplet  transitions  using  the  TDDFT  methodology. 
The  LDA  and  GGA  calculations  can  be  efficiently 
performed  by  using  density  fitting  thus  replacing 
formally  N4  integral  algorithm  with  the  more  efficient  N3 
method.  We  found  that  the  electronic  transitions  did  not 
depend  significantly  on  the  particular  DFT  functional 
approximation  and  the  B3LYP  method  was  therefore 
used  in  the  majority  of  our  calculations.  We  have  also 
used  the  configuration  interaction  method  with  single 
excitations  (CIS)  within  the  semiempirical 
ZINDO/IND02  approach. 

3.  RESULTS 

A  level  of  methodology,  which  will  be  chemically 
meaningful  to  guide  experiment  and  push  forward  our 
understanding  of  these  systems,  was  determined  by 
evaluating  structural  and  property  invariance  with 
increasing  sophistication  of  basis  set  and  treatment  of 


dynamic  electron  correlation.  Thus  sensitivity  of  results 
to  the  choice  of  quantum  chemistry  model  will  be 
discussed  first.  Next,  results  for  the  PC  with  the  PEG 
attached  will  be  presented  and  compared  with  the  recent 
experimental  findings  (Huang  et  al,  2004).  The  effect  of 
stacking  PC  rings  is  then  explored  in  section  3.3. 

3.1  Tests  of  quantum  model 

Attaching  two  PEG  units  with  a  mass  of  75  to  the 
SiPC  framework  allows  us  to  test  the  performance  of 
various  DFT  functionals  on  structural  and  property 
predictions.  The  structure  of  SiPC-PEG75  (Fig.  1)  was 
optimized  in  the  gas  phase,  followed  by  UV-Vis  and 
RSA  property  analysis,  as  described  in  section  2.  Table  1 
shows  the  main  excitation  energies  (in  nm)  and  oscillator 
strengths  (f)  calculated  at  the  following  levels  of  theory: 
B3LYP/631G*,  B3LYP/631 1G*,  LDA/631G*, 

PW91/631G*,  PBE/631G*,  PBE/6311G*,  including  also 
results  for  the  ZIND0/IND02  approximation. 


Figure  1.  Structure  of  SiPC-PEG75 


The  size  of  the  ZINDO  CIS  expansion  is  given  in 
parentheses  in  Table  1.  In  addition  to  structural  results 
performed  at  the  PBE/DNP  level  of  theory,  ZINDO 
calculations  were  also  performed  using  semiempirical 


Table  1.  Main  excitation  energies  (AE  in  nm)  and 
oscillator  strengths  (f)  from  UV-Vis  and  RSA 
calculations  of  the  SiPC-PEG75  model  complex. 

Method\Excited  state  Q  B  R 

DFT 


B3LYP/631G* 

613(0.37) 

335(0.54) 

479(0.31) 

B3LYP/6311G* 

618(0.38) 

338(0.61) 

487(0.30) 

LDA/631G* 

624(0.26) 

368(0.21) 

585(0.41) 

PW91/631G* 

626(0.27) 

371(0.20) 

578(0.39) 

PBE/631G* 

625(0.27) 

371(0.24) 

579(0.30) 

PBE/6311G* 

628(0.28) 

362(0.23) 

582(0.41) 

ZINDO/IND02 

ZINDO(5)//DFT 

671(0.78) 

313(2.40) 

ZINDO(  1 0)//DFT 

671(0.78) 

317(2.07) 

ZINDO(20)//DFT 

698(0.77) 

337(1.33) 

ZINDO(3  0)//DFT 

705(0.74) 

340(1.27) 

ZINDO(40)//DFT 

717(0.69) 

343(1.28) 

ZINDO(50)//DFT 

728(0.66) 

349(1.23) 

ZINDO(  1 0)//PM3 

658(0.72) 

354(1.07) 

ZINDO(  1 0)//AM  1 

643(0.68) 

367(0.99) 

PM3  and  AMI  optimized  structure  of  the  complex.  The 
excitations  with  the  largest  value  of  oscillator  strength,  f 
from  Q  and  B  bands  are  presented  together  with  the  main 
RSA  transition  (R).  Currently,  the  ZINDO  program  does 
not  have  the  capability  to  calculate  RSA  spectra.  The 
excitation  energies  of  local  and  gradient-corrected  DFT 
approach  appear  to  be  red  shifted  in  comparison  to  the 
adiabatic,  B3FYP  method.  The  R  band,  in  particular  is 
red  shifted  by  about  60-100  nm.  Effects  of  basis  set 
extension,  including  diffuse  functionality  is  not 
significant  with  a  shift  in  the  spectrum  by  about  10  nm 
with  the  more  extensive  basis  set.  We  have  found  the 
density  fitting  algorithm  to  be  computationally  much 
more  efficient  while  not  compromising  the  quality  of  the 
results  for  electronic  spectra.  Therefore  results  obtained 
at  the  EDA  and  GGA  level  of  DFT  approximation  were 
performed  using  the  standard  density  fitting  algorithm  as 
implemented  in  the  GAUSSIAN03  program. 

The  character  of  transitions  was  studied  (see  Table 
2)  and  it  was  found  to  follow  approximately  the 
Goutermann’s  assignments.  The  Q  band  is  mainly  due  to 


Table  2.  Character  of  transitions  for  major  excitation 
energies  from  B3LYP  and  PBE  calculations  with 
631 1G*  basis  for  SiPC-PEG75  model  complex. 

All  electronic  states  belong  to  !A  irrep. 


AE  f 


transition 


B3LYP: 

618  0.375 

605  0.371 

339  0.092 

338  0.609 

337  0.524 

321  0.120 

320  0.270 

622  0.141 

528  0.241 

516  0.115 

487  0.300 

483  0.113 


Singlet— ►singlet 
H— ►L 
H— ►L+l 
H-6— ►L 
H-6— ►L 
H-6— ►L+l 
H— >L+5 
H-14— ►L+l 
Triplet— ►Triplet 

H-l— ►L(A) 
H-8— ►L(B) 
H— >L+6(A) 
H-9— ►L(B) 
H-13— ►L(B) 


weight 


0.603 

0.606 

0.595 

0.577 

0.489 

0.509 

0.501 

0.586 

0.882 

0.696 

0.514 

0.740 


PBE: 

Singlet- 

-►singlet 

648 

0.195 

H— ►L 

628 

0.284 

H— ►L+l 

615 

0.093 

H-3— ►L 

404 

0.119 

H-13— ►L 

372 

0.122 

H-17— ►L+l 

371 

0.119 

H-17— ►L+l 

363 

0.141 

H— >L+6 

362 

0.224 

H— >L+4 

Triplet- 

-►Triplet 

674 

0.111 

H-9— ►L(B) 

582 

0.412 

H-l  1— ►L(B) 

0.483 

0.578 

0.616 

0.495 

0.383 

0.471 

0.539 

0.406 

0.716 

0.467 


the  HOMO(H)  to  LUMO  (L)  and  L+l  excitations.  The  B 


a) 

Figure  2.  a)  Unit  cell  of  SiPC-PEG25o 


b)  c) 

with  500  water  molecules.  Conformers  of  b)  SiPC-PEG75  and  c)  SiPC-PEG25o 


band  excitations  can  be  interpreted  as  excitations  from 
several  levels  below  H  to  L  and  L+l  levels.  The  B  band 
appears  to  be  more  sensitive  to  the  level  of  DFT 
approximation  than  the  B  band.  The  PBE-originated  Q 
band  is  broader  and  at  the  maximum  corresponds  to 
transitions  to  higher  virtual  orbitals  than  in  the  case  of 
the  hybrid  B3LYP  functional  (see  Table  2).  The  RSA 
spectra  differ  significantly  depending  on  selection  of  the 
DFTfunctional.  The  B3LYP  band  is  broad  and  covers  the 
entire  green  part  of  the  spectrum,  while  the  PBE  RSA 
band  is  very  narrow  and  shifted  to  the  red.  The  character 
of  the  transition  is  also  different  with  alpha  electrons 
(from  the  spin  polarized  calculations)  participating  in  the 
spectra  for  the  B3LYP,  while  the  beta  electrons  dominate 
the  PBE  transitions. 

Both  the  local  PBE  functional  as  well  as  the  B3LYP 
hybrid  functional  approaches  predict  the  absorption  of 
the  red  light  and  hence  a  complementary,  green  color  of 
the  complex.  The  B3LYP  method  shows  the  strong  RSA 
absorption  in  the  expected  green  region,  while  the  GGA- 
type  method  shows  triplet-triplet  transitions  only  in  the 
orange  part  of  the  spectrum.  Based  on  these  results  we 
use  the  PBE  method  to  determine  the  color  of  the 
complex,  while  the  B3LYP  method  is  used  to  determine 
the  RSA  properties  of  the  system.  As  such  the 
determination  of  the  UV-Vis  spectra  can  be  done  with  a 
computationally  much  more  efficient,  PBE,  method. 

The  semiempirical  ZINDO  method  offers  extremely 
fast  calculations  of  UV-Vis  spectra.  Both  the  Q  and  B 
band  become  red  shifted  with  an  increase  in  the  CIS 
ZINDO  expansion.  The  AMI  and  PM3  -optimized 
structures  differ  significantly  from  the  DFT  structure  and 
therefore  are  not  recommended  for  study  of  PC 
complexes.  The  ZIND0/IND02  parameters  at  the  DFT 
geometry,  however,  reproduces  exceptionally  well  the 
experimental  data  (Huang  et  al.,  2004)  at  low  values  of 
CIS  expansion,  thus  making  it  an  attractive  method  for 
fast  evaluation  of  the  UV-Vis  spectra. 


Following  the  procedure  described  in  Section  2, 
several  structures  of  SiPC-PEG  structures  with  PEG 
mass  of  75  and  250  were  generated  using  the  MD 
approach  in  the  presence  of  water.  Next,  structures  were 
optimized  and  the  energy  conformers  differing 
significantly  in  energetic  values  were  used  for 
subsequent  excited  states  calculations.  Fig.  2a,  presents  a 
snapshot  from  the  MD  computation  and  Fig.  2b 
illustrates  the  SiPC-PEG75  conformer  that  differs  in 
energy  by  2.0  kcal/mol  from  that  presented  in  Fig.  1.  Fig. 
2c  presents  two  SiPC-PEG250  conformers  differing  in 
energy  of  4.1  kcal/mol.  In  addition,  snapshots  from  MD 
simulations  were  selected  with  several  (3)  water 
molecules  within  5  A  of  the  SiPC-PEG75  conformer. 

Table  3,  presents  the  main  excitation  energies  for 
selected  conformers,  cl  and  c2.  The  B3LYP/6- 
31G*//PBE/DNP  level  of  theory  was  used.  The 
ZINDO(CIS=  1 0)/INDO2//PBE/DNP  results  are  also 
reported  in  Table  3  for  selected  conformers.  Comparison 
of  major  excitation  peaks  indicates  that  the  conformation 
of  PEG  chain  does  not  affect  either  the  UV-Vis  or  RSA 
spectra.  Including  explicit  water  molecules  also  does  not 


Table  3.  Main  excitation  energies  (AE  in  nm)  and 
oscillator  strength  (f)  from  UV-Vis  and  RSA 
calculations  for  SiPC-PEG  (X=75,250)  complexes. 


System\Excited  state 

Q 

B 

R 

DFT 

cl  SiPC-PEG75 

613(0.37) 

335(0.54) 

479(0.31) 

c2  SiPC-PEG75 

613(0.36) 

337(0.68) 

477(0.22) 

cl  SiPC-PEG75(3w) 

611(0.36) 

334(0.45) 

516(0.20) 

cl  SiPC-PEG25o 

611(0.35) 

337(0.47) 

485(0.25) 

c2  SiPC-PEG25o 

ZINDO(CIS= 

611(0.34) 

10)/INDO2 

338(0.42) 

500(0.18) 

c2  SiPC-PEG75 

671(0.78) 

327(0.56) 

cl  SiPC-PEG25o 

671(0.78) 

315(1.66) 

c2  SiPC-PEG250 

670(0.78) 

316(1.46) 

3.1  Excited  states  of  SiPC-PEGx(x=75, 250)  appear  to  significantly  affect  the  spectral  results.  The 

spectra  obtained  with  the  short  PEG  chain  (75)  and  the 
long  PEG  chain  (250)  are  also  very  similar.  This  is  in 


a)  b) 

Figure  3.  Conformers  of  SiPC  with  PEG75  and  Cl. 


c) 


line  with  experimental  measurements,  obtained  for  a 
longer  PEG  chains  with  mass  of  550  and  750  (Hunag  et 
al.,  2004).  Experimental  results  indicate  an  intense  and 
sharp  Q  band  at  about  671  nm  and  a  B  band  at  about  354 
nm.  These  values  can  be  compared  to  our  findings  of 
about  620  and  340  nm,  respectively.  Although  this  level 
of  DFT  calculations  leads  to  a  blue  shifted  Q  band,  the  B 
band  matches  closely  the  experimental  results.  The 
ZINDO  spectrum  matches  closely  to  the  experimental 
data,  particularly  for  the  Q  band. 

We  conclude  that  the  UV-Vis  and  RSA  spectral 
results  are  relatively  insensitive  to  influences  of  the 
environment,  as  well  as  to  the  size  and  specific 
conformation  of  PEG.  We  confirm  that  the  SiPC-PEG 
complex  is  green  in  color,  and  has  primary  RSA 
absorption  in  the  visible  range  of  the  UV-Vis  spectrum. 

3.2  Effect  of  SiPC  stacking  on  spectra 

To  elucidate  further  the  impact  of  PEG  on  spectral 
properties,  we  performed  calculations  for  various 
conformers  of  SiPC  without  PEG  attached,  with  PEG 
attached  on  one  side  only,  and  including  other 
substituents  such  as  chlorine.  The  primary  motivation  for 
this  investigation  was  to  assess  a  possible  color  changes 
due  to  incomplete  attachment  of  PEG  fragments  that  may 
result  in  stacking  of  several  PC  rings.  Fig.  3  illustrates 
several  structures  with  incomplete  PEG  attachment, 
while  Figure  4  (b-d)  shows  optimized  structures  used  to 
study  effects  of  stacking.  Table  4  lists  UV-Vis  spectral 
results  in  the  visual  range  of  400  to  700  nm,  which 


includes  the  Q  band  of  the  PC  complexes.  The  B3LYP/6- 
31G*//PBE/DNP  DFT  level  of  theory  was  used  and  the 
ZINDO(CIS=30)/INDO2//PBE/DNP  semiempirical 
method  was  utilized  for  comparison. 

The  trends  in  absorption  spectrum  as  displayed  in 
Table  4  are  similar  for  both  the  DFT  and  ZINDO  results, 
although  the  absolute  values  of  excitation  energies  are 
vastly  different.  Apparently  substituting  PEG  oligomers 
with  chlorine  atom  (Fig.  3c,  Fig.  3d  and  Fig.  2b)  does  not 
affect  significantly  the  position  of  the  Q  band. 
Additionally,  water  molecules  forming  hydrogen  bonds 
with  the  complex  (structure  F4a  from  Fig.  4a)  also  does 
not  appear  to  significantly  alter  the  spectral  properties. 


Table  4.  Main  excitation  energies  (AE  in  nm)  and 
oscillator  strength  (f)  in  the  visual  range  as  calculated 
from  DFT  and  ZINDOmethods.  Structures  are  shown 
in  Fig.3  (F3)  and  Fig.4  (F4).  Detailed  explanation 
given  in  the  text. 

System\Excited  state 


DFT 

ZINDO 

F3a(D4h) 

516(0.67) 

450(2.69) 

F3a(fromF3b) 

667(0.29) 

780(0.57) 

F3b 

658(0.31) 

774(0.66) 

F3c 

620(0.36) 

717(0.73) 

F3d 

617(0.36) 

713(0.73) 

F4a 

626(0.36) 

710(0.74) 

F4b 

523(0.77) 

517(1.82) 

F4c 

536(0.60) 

497(1.71) 

Figure  4.  Structure  of  SiPC-PEG75  with  hydrogen-bonded  water  molecules  and  complexes  of  stacked  SiPC  fragments. 


Distortion  of  the  PC  ring,  such  as  displayed  in  Fig 
3b,  shifts  the  spectrum  towards  the  red.  Structure  F3b 
without  the  PEG  oligomers  (denoted  as  F3a  (fromF3b )  in 
Table  4)  does  not  change  the  spectral  results  of  the 
original  F3b  structure  indicating  that  the  PEG  oligomers 
do  not  influence  significantly  the  Q  band.  However,  the 
Si-PC  complex  optimized  under  the  D4h  symmetry 
constraints  shows  a  marked  difference  in  spectral  results. 
At  the  B3LYP/6-31G(d)  level  of  theory,  one  finds 
several  absorption  energies  in  the  green  region,  with  a 
peak  at  516  nm.  The  high  symmetry  of  the  SiPC 
structure  leads  to  a  much  shorter  Si-N  bond,  by  almost 
0.1  A,  as  compared  to  structures  with  the  PEG  attached. 
The  net  charge  on  the  central  Si  atom,  according  to  the 
Hirschfeld  analysis,  is  also  markedly  different  and 
amounts  to  0.55  as  compared  to  about  0.42  for  the  F3b 
structure  with  PEG  oligomers. 

The  effect  of  stacking  has  a  significant  effect  on  the 
color  of  the  molecular  system.  Results  in  Table  4  show  a 
significant  blue  shift  in  the  spectrum  for  structures  F4b, 
F4c  and  F4d  (see  Fig.  4).  We  note  that  the  structure 
represented  in  F4d  has  yet  to  be  calculated  due  to  its  size. 
Without  the  PEG  oligomers  attached,  the  SiPC  rings  tend 
to  form  an  internal  Si-Si  bond,  while  for  all  stacked 
structures  a  strong  n-n  interaction  is  expected. 

We  conclude  that  the  PEG  oligomers  affect  the 
spectroscopy  of  SiPC  system  by  imposing  geometrical 
constraints  and  by  minimizing  stacking  in  the  SiPC 
complexes.  Upon  attachment  to  the  SiPC  ring,  PEG 
causes  local  deformation  of  the  ring,  increasing  the  Si-N 
distance  and,  in  the  case  of  non-symmetrical  complexes 
such  as  F3d,  extending  the  Si  out  of  the  PC  plane.  In  the 
absence  of  axial  substituents  the  Si  atom  transfers  charge 
to  the  neighboring  nitrogen  atoms.  However,  since  both 
PEG  and  Cl  atom  withdraw  electrons  from  the  SiPC 
system  and  the  net  charge  on  the  Si  atom  decreases,  it 
indicates  the  charge  flows  from  the  PC  ring  through  Si  to 
the  ligands.  This  rearrangement  of  electronic  structure 
because  of  underlying  geometrical  changes  causes  the 
red  shift  of  the  Q  band,  thus  making  SiPC-PEG  material 
transparent  in  the  green  part  of  the  visual  range. 


CONCLUSIONS 

In  this  work,  we  have  evaluated  various  quantum 
chemistry  methods  for  calculations  involving  structure 
and  spectral  results  of  a  series  of  PEG  complexes.  In 
particular,  UV-Vis  and  RSA  spectra  of  Si- 
phthalocyanine  complexes  with  axially  attached  PEG 
poly(ethylene  oxide),  have  been  investigated.  Based  on 
these  computational  studies,  we  find  that  GGA  and  LDA 
density  functional  approaches,  using  a  computationally 
efficient  algorithm  of  fitting  density,  reproduces  well  the 
UV-Vis  spectrum  found  experimentally  for  complexes  of 


this  type.  The  ZINDO  semiempirical  method  predicts 
very  accurate  UV-Vis  spectra  provided  that  the  CIS 
expansion  is  short,  and  the  underlying  optimized 
structure  is  optimized  using  higher  order  quantum 
chemical  methods,  such  as  provided  by  density 
functional  theory.  Reliable  results  for  RSA  spectra  can 
only  be  obtained  using  higher  level  computational 
methods,  such  as  provided  by  the  hybrid  B3LYP  method. 
The  DFT  methods  require  a  modest  basis  set  expansion, 
however,  further  study,  including  diffuse  functions  is 
ongoing  to  confirm  this  assessment.  These  conclusions 
are  important  for  the  study  of  much  larger  complexes 
including  napthalocyanines  and  their  derivatives. 

In  this  study,  we  have  found  that  the  length  of  PEG 
oligomers,  the  particular  conformational  preference,  as 
well  as  the  water  environment  simulated  with  inclusion 
of  explicit  water  molecules,  does  not  significantly  affect 
either  the  UV-Vis  or  RSA  spectra.  However,  partial 
attachment  of  PEG  oligomers  to  SiPC  complex  may  lead 
to  stacking  of  SiPC  rings,  significantly  affecting  UV-Vis 
spectral  properties.  The  system  will  thereby  no  longer  be 
transparent  in  the  desirable  green  range  of  the  spectrum. 

The  results  of  this  paper  improve  our  understanding 
of  the  binding  of  PEG  oligomers  to  SiPC  complex  and 
also  may  help  in  selecting  axial  substituents  for  synthesis 
of  phthalocyanines  with  desirable  optical  properties. 
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